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Abstract 

We show that scattering amplitudes between initial wave packet states and cer- 
tain coherent final states can be computed in a systematic weak coupling ex- 
pansion about classical solutions satisfying initial value conditions. The initial 
value conditions are such as to make the solution of the classical field equations 
amenable to numerical methods. We propose a practical procedure for comput- 
ing classical solutions which contribute to high energy two particle scattering 
amplitudes. We consider in this regard the implications of a recent numerical 
simulation in classical SU(2) Yang-Mills theory for multiparticle scattering in 
quantum gauge theories and speculate on its generalization to electroweak theory. 
The generalization of our results to complex trajectories allows its application to 
any wave packet to coherent state transition. Finally, we discuss the relevance of 
these results to the issues of baryon number violation and multiparticle scattering 
at high energies. 
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1 Introduction 



The behavior of quantum scattering amplitudes at high center of mass energies, E ~ E* = 
m/g 2 , is an open question currently under investigation. Here, m is a mass scale charac- 
terizing the particle excitations of the theory, and g is a small coupling constant controlling 
the nonlinear interactions. It has been speculated that multiparticle scattering from initial 
states with only a small number of particles could be unsuppressed at these energies. A 
second related issue is whether anomalous processes, such as baryon number violation in 
electroweak theory where E* ~ O (M w /a w ) ~ 10 TeV, are unsuppressed at energies of the 
same order JTJ. 

Standard perturbation theory is clearly inadequate to address either of these issues. Re- 
cent investigations have focused on expansions about nonperturbative solutions of classical 
field equations, such as instantons [0], ||. However, none of the previously proposed solu- 
tions provides a stationary point for a controlled semiclassical expansion at the energies of 
interest, E ~ E* ||. In general, the semiclassical expansion about the proposed instanton- 
like saddlepoints is controlled only at low energies and provides little or no information at 
E ~ E*. 

In this paper, we propose a method for finding classical solutions which will prove suitable 
saddlepoints of scattering amplitudes even at E ~ E*. The equations and boundary condi- 
tions which must be satisfied are dramatically simplified compared to a previous approach to 
the same problem by Mattis, McLerran and Yaffe [[|. In particular, the resulting equations 
and boundary conditions for the classical saddlepoint are no longer of the integro-differential 
type. Rather, the equations of motion are source-free, 

5S 
6<f>(x)~ 

while the boundary conditions are of the initial value type, which can be integrated forward 
in time using standard numerical algorithms. This is also a departure from the approach of 
Khlebnikov, Rubakov and Tinyakov ^j, who derive similar source-free equations, but for 
an initial coherent state, rather than a two particle state. 

We discuss the implications of some recent numerical simulations in classical field theory 
for high energy two-particle scattering. In particular, in 4-dimensional pure SU(2) Yang- 
Mills theory, Gong, Matinyan, Miiller and Trayanov || have found an instability of high 
frequency standing waves of pure glue to decay to low frequency modes. When interpreted 
within our formalism, this result leads to the speculation that multiparticle scattering in 
SU{2) gauge theories may be unsuppressed (at least exponentially) at high energies. We 
give a generalization of our method to complex trajectories which allows their application 
to essentially any wave packet - coherent state scattering amplitude. We finally discuss how 
our formulation relates to some previous work by other authors on high energy multiparticle 
and anomalous processes in electroweak theory || f| . 
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2 Formalism 



We wish to calculate, in the stationary phase approximation, the on-shell S-matrix element 
between an initial two particle state and an arbitrary final state 0. The result will be a 
boundary value problem which determines the classical field which is the stationary point of 
the S-matrix element. 

To this end, we express the kernel of the S-matrix in a basis of coherent states, first used 
in this context in Q. The initial and final states are defined by sets of complex variables 
a = {a k }, b* = {& k }, respectively. A coherent state |a k ) is an eigenstate of the annihilation 
operator d k : d k |a k ) = a kl a k)- The relevance of coherent states for the scattering problem 
stems from the fact that upon differentiation they generate all momentum eigenstates (see 
eq. j2.lUI below). 

For simplicity, we illustrate the case of a single real massive scalar field with self-coupling 
g 2 in 3 + 1 dimensions. The results can be generalized to the case of gauge fields. First, 
we express the transition amplitude from an initial coherent state |a) at time Tj to a final 
coherent state |b*) at time Tf, in terms of the path integral by inserting a complete set of 
eigenstates of the field operator on the initial and final time slices (see figure 1) 

(b*|C/|a) = Jd^id<f> f (b*|0/) (<t> f \U\<l) i ) <&|a) . (2.1) 

U is the evolution operator between time Tj and Tf. A "position" eigenstate of the field 
operator <ft is denoted \<p) and <pij = 0(Tj). Then, from ( |2.1| ), we obtain the S-matrix 
kernel in a compact form in terms of path integrals 

(b*|£|a) = S[b*,a] = lim / c% dfc e B f e Bl [' e lSW , (2.2) 

Ti,Tf -*^oo J J 

where S [<p\ is the action functional. The path integral appearing here is over fields obeying 
the boundary conditions 0(Tj) = 4>ij- The functional Bf is 

B f [b*,<f>f} = -y d *kb^ k e 2 ^ T f-yd 3 k^<P f (k)M-k) (2.3) 
+ J d 3 k V2~J^ e iuJkT f b* k 4> f {-k) , 

in terms of which the wave functional of the final coherent state is 

(b*|^) = exp( J B / [b*,0 / ]) . (2.4) 

Similarly, the functional Bi is 

Bi [a, ^ = -~J d 3 k a k a_ k e" 2 ^ Tl -\jd 3 k ^(^^(-fe) (2.5) 
+ J d 3 k e~ M a k ^(k) 

1 The authors of || extremize momentum-space Greens functions. See Section 6. 
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is related to the wave functional of the initial coherent state 

( (j>i | a ) = exp ( Bi [a, ) . 

Here the 3-dimensional Fourier transform is defined 

d 3 x 



Ak-x 



(2vr) 3 / 2 



(2.6) 



(2.7) 



and is related to the residue of the 4- dimensional Fourier transform when the field reduces 
to a plane wave superposition as | Tjj | — > oo 0. 

t 



Figure 1: Spacetime, with two slices. 



The kernel (|2.2|) is a generating functional for S-matrix elements between any initial and 
final N particle states, by functional differentiation with respect to arbitrary and We 
now use this fact to construct a kernel for scattering from initial two particle states. We 
define an initial two particle (wave packet) state at t = Tj 



p,-p) 



d 3 k a R (k) al / d 3 k' a L (k') at | 0) 



(2.8) 



x — * — * 

where is a creation operator, and ctR^(k) are arbitrary smearing functions of k, localized 
around some reference momenta p and — p respectively. The wave packets are normalized so 
that! 

(2.9) 



J d 3 k | a RtL (k) 



The regime relevent to the high energy multiparticle scattering problem is where, in the 
limit g 2 — > 0, the wave packet is localized around momenta \p\ ~ m/g 2 with a characteristic 
width, Ap ~ m. 



2 We use the non-relativistically normalized commutator 



a k , , 



S 3 (k - £') . 
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This state can be generated by functional differentiation of the coherent state |a) with 
respect to Ok 



\p, -p) = / d 3 k d 3 k' a R (k) a L (k') - — - — | a) 
J oa\c oa\ci 



(2.10) 

a = 



So, differentiating under the functional integral, the S-matrix element between the two 
particle state ( |2.10| ) and any final state involves the following functional at t = Ti 

5 6 



exp ( Bi [a, 



a=0 



(2.11) 



2 v /^ k 7 e -^ k +^m exp ( Bi [0> 



after dropping a term which vanishes in the limit Ti — > — oo. The last factor here is simply 
the normalization of the initial position eigenstate 

exp (~j d 3 k u k 4>i{k) ) . (2.12) 

We combine this with the smearing functions and finally obtain an S-matrix kernel for 
the scattering of two wave packets into arbitrary final states, 

S [b*, 2] = lim / (fyfdfc a R -<pia L - & e B ^ b M + m[0 ^' ] D(j) e i5[</,] , (2.13) 

where we have denoted the initial state (|2.10|) by "2". Here we have used the following 
compact notation for the initial state factors 



a 



d 3 k a{k) <Pi{k) e~' M . (2.14) 



We now derive the boundary value problem obeyed by the classical field in the stationary 
phase approximation of (|2.13| ). The stationary phase approximation of ordinary integrals 
suggests that the initial state factors be included in the extremization of the path integral i. 
These factors will certainly impact the stationary phase in the kinematic regime relevant 
to the high energy scattering problem, where the initial wave packets are peaked around 
\p\ ~ m/g 2 , as expansions around the instanton indicate More generally though, we will 
show in Section 4 that this procedure provides a consistent weak coupling expansion at any 
fixed center of mass energy. 

We include initial state factors in a straightforward manner by first exponentiating the 
initial state factors into an "effective action" , so that 

S[b*,2) = lim / dfodfaDcj) e r , (2.15) 



3 The Stirling approximation to the Gamma function is a classic example. 
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where the effective action T is 



T [0] = lna R -faa L -fa + B t [0,fa] + iS[4>] + B f [b*,(f> f } , (2.16) 

after dropping a term which vanishes as Tj — > — oo. 

We can now derive the boundary value problem by varying the effective action. Varying 
the entire exponent T with respect to <p(x) for Tj < t < Tf gives the source-free equations of 
motion 

m = (2 - 17) 

Varying the entire exponent with respect to fa(k), gives 



ifa(k) + u k fa(k) = V2^ ( + 



ot R ■ (pi a L - 



(2.18) 



The first term on the left hand side comes from a surface term in the action S. The other 
terms come from variation of the wave functional at £ = Tj. This boundary condition involves 
both the positive and negative frequency parts of the field, unlike the boundary condition 
which arises for an initial coherent state |a), which depends on the negative frequency com- 
ponent of 4> only f2j. 



The boundary condition ( 2.18 ) at the initial time slice is rather complicated. However, it 
can be simplified since a real field <p may be written in the asymptotic region t — Tj — ► — oo 
as a plane wave superposition 

fa(k) = -=L (« k e-^ + «* k e^ Ti ) . (2.19) 
Equation ( |2.18| ) then reduces to the requirement 



a R (k) + a L (k) 

"k = t^T72 ' ( 2 " 2 °) 

+ J d 3 k aR^k) aL{k) 



using the normalization ( |2.9| ). This solution is consistent with physical intuition, the classical 
field reducing to the initial particles at early times. The overlap of the left- and right- moving 
wave packets in the denominator is very small for narrow high energy wave packets. 

If the field is real, then the negative frequency part equals the complex conjugate of the 
positive frequency part, and the field fa is determined as in ( |2.19| ) and ( |2.20| ). The real 
initial condition can be integrated forward, and uniquely determines the final field 0/. We 
say more about the case of complex stationary points of real fields in Section 5. 

It remains to be seen in what sense it is a good approximation to include the initial 
state prefactors in the stationary phase calculation. We will demonstrate in Section 4 that 
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the stationary phase solution which results from this approach provides a controlled weak 
coupling expansion of the scattering amplitude. 
Finally, varying with respect to <pf{k) gives 

- i<j>f{k) + co k ^ f (k) = v/2^6I k e^ T / . (2.21) 

Again with a free field form as in ( [2.19| ), this boundary condition may be re-expressed as 

<P f (k) = -=L= (b k e~^ T f + b*_ k e^ T f) . (2.22) 

We see that the precise scattering amplitude which is extremized is only determined at the 
end of the calculation, by the asymptotic form of the classical solution in the far future, 
Tf — > oo. For the initial conditions given by ( p. 18 ), the scattering amplitude corresponds 



to a transition from wave packet states to a final coherent state described by \b*) from 
( |2.21| ). Note that neither a wave packet state @ nor a coherent state are eigenstates of the 
Hamiltonian, so that neither have a definite energy. 

When the classical field satisfies the equation of motion ( |2.17| ) and the boundary condi- 
tions (p. 18 ) and ( 2.21 ), there are no initial and final state corrections to leading order in the 



semiclassical expansion. This is not the case for the instanton, which satisfies the equation 
of motion but not the correct boundary conditions. Instead, there are linear terms in the 
fluctuation about the instanton and these generate initial and final state corrections. We 
will discuss this further in Sections 4. 



3 The Program 

The results of the previous section establish an important fact: 

The classical field obtained by starting with two incident wave packets and evolv- 
ing them forward using ( |2.17| ) is the dominant contribution to some two particle 
scattering amplitude, in the semiclassical approximation. 

In the case of a purely Minkowskian classical solution, the result for the corresponding scat- 
tering amplitude will not contain an exponential suppression. Thus, if a real time classical 
trajectory can be found which connects wave packet initial conditions to an "interesting" 
final state, for example a multiparticle or fermion number violating state, our results imply 
that the corresponding scattering amplitude is unsuppressed. Classical solutions which ei- 
ther exist in complex time or are themselves complex are relevant to classically forbidden 
transitions and will be discussed in Section 5. 

4 Except in the limit of au(k) ~ S 3 (k — p), where it reduces to a plane wave. 
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The boundary conditions given by (|2.19|) and Q2.20| ) are amenable to straightforward 
numerical integration. The boundary conditions specify <p(x, TA and 0(x, T;), which are 
sufficient to construct (p(x, t) for all subsequent t > Ti, given a discretization of the equations 
of motion fl2.17p. In principle, it is possible to generate an infinite number of classical 



trajectories, each relevant to a particular initial wave packet scattering amplitude. 

Of course, it is not clear a priori that the final state which results from the classical evolu- 
tion will be one which is interesting to the problem of multiparticle production or anomalous 
baryon number violation. However, the above observation summarizes the relevance for 
scattering amplitudes of computations by Rajagopal and Turok |J, and also Goldberg, Nash 
and Vaughn f7j. Rajagopal and Turok studied the classical scattering of wave packets in the 
Abelian Higgs model, while Goldberg et al. studied the classical dynamics of 4 theory. Nei- 
ther group found that energy was readily transferred from high frequency to low frequency 
modes, a signal which would indicate the production of a final state with many (soft) parti- 
cles. Therefore, the scattering amplitudes which are dominated by their classical trajectories 
are probably not relevant to multiparticle production or fermion number violation. 

On the other hand, the initial indications on the behavior of non-Abelian classical tra- 
jectories seem more promising. Some recent results on the classical behavior of pure SU(2) 
Yang-Mills theory by Gong, Matinyan, Miiller and Trayanov || may provide insight into so- 
lutions of the boundary value problem presented above. These authors have considered the 
classical stability of a stationary mono-color wave in Yang- Mills theory (in A$ = gauge): 

A^(x,t) = 5 i3 5 c3 A cosk x cosu t, (3.1) 

where c is a color and i is a spatial index i. Small amplitude variations of the field in 
directions of different color are found to lead to an instability with long wavelength. 

The Hamiltonian evolution of a perturbed standing wave is equivalent to the evolution 
of the initial conditions relevant to the scattering of two perturbed plane waves in SU(2) 
gauge theory. The discovery of an instability of the initial configuration to decay into long 
wavelength modes implies the existence of a classical trajectory connecting initial high energy 
plane waves to long wavelength modes in the final state. It therefore suggests that the 
corresponding 2 — > many gluon scattering amplitude may be unsuppressed ! Unfortunately, 
the relevance of weak coupling calculations to pure gauge theory is unclear, due to the 
asymptotic freedom of the theory. It is possible that arguments similar to those from deep 
inelastic scattering may be used to justify the semiclassical result. In the limit that all of the 
energy scales (including the energies of the many "soft" outgoing gluons) are large compared 
to the intrinsic mass scale of the theory Asu(2), it seems plausible that the relevant running 
coupling constant is small. It is also important to determine whether the instability persists 
for wave packet initial conditions, as the plane wave limit is never achieved at an actual 
accelerator. 



5 NB - A travelling wave is stable in Yang-Mills theory 
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Another class of trajectories in pure SU(2) Yang-Mills theory are those found by Farhi 
et al. [[|. These solutions in the spherical ansatz correspond to spherical shells of glue 
that originate at infinity, collapse inward, and bounce back to infinity, leaving behind some 
fractional topological charge. There are some open questions concerning such solutions - in 
particular, their stability to perturbations which are outside the spherical ansatz, and their 
relevance to fermion number violation. However, in the absence of fermions, they can be 
interpreted within our formalism as extremizing scattering amplitudes between spherically 
symmetric initial and final quantum states. 

Of course, our eventual goal is to address the electroweak theory, where the Higgs mecha- 
nism cuts off the infrared growth of the coupling constant. It would be extremely interesting 
if simulations of wave packet scattering could be performed in a spontaneously broken SU(2) 
gauge theory. 

Clearly, there is much numerical work to be done to address these issues. 



4 Small Fluctuations 

In this section, we construct the perturbative expansion around the classical solution of the 
BVP. We will demonstrate that corrections to the leading contribution are systematically 
suppressed in powers of the small coupling constant for any fixed energy. The resulting small 
coupling expansion differs from the usual one (e.g.- instanton) and the subtleties involved 
will be discussed at the end of this section. 

As usual, we proceed by expanding the quantum field around a classical background 
field C , which is a solution to the classical boundary value problem derived in Section 2. 

<j){x) = (f) c (x;/j) + u{x) . (4.1) 

The classical solution depends in general on a collection of collective coordinates denoted 
here by /i, corresponding to the invariances of the scattering amplitude || pT0| . We expect 
invariances of translation and scale size to be broken by the choice of cxR^{k)] however, fi 
may include gauge coordinates in a gauge theory. Equation (|4.1| ) involves an expansion on 
the initial and final time slices, which we can express in momentum space as 

= 0) + v ij$) for t = T ui ( 4 - 2 ) 



where $ j(k) is known explicitly for any given initial and final wave packets, from ( |2.19| ) and 



Note that the classical solution is of order O(g ), since it is non-trivial in the limit 
as g — > 0. In fact, when g vanishes, the background field is simply the initial condition 
propagated forward in time with the free Hamiltonian, 

c (t,fc) = (lite^ 1 + ul k e-**) , (4.3) 
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where «k is as in ( |2.20|) . In this limit, the two wave packets simply pass through each 
other. For non-vanishing g, the classical field <p c will have some complicated dependence on 
g, reflecting the non-linearity of the theory and the external boundary conditions. 

The kernel of the S-matrix ( |2.13| ) is composed of three functional integrals which we write 

as 



S[b*,2) 



dfa dfa a R .<j>ia L - <\> t e B ^^ + B ^ 



iS[<f>] 



(4.4) 



where the limit \Tij\ — > oo is understood. Now we formally expand each of the terms in 
( [4.4p using ( f4.1| ). The action becomes 



S^ + v] = S [</> c ] + I d'x v(x) ^(x) + S 2 [v] + S int [u] 



(4.5) 



after integrating by parts and retaining the surface terms in the time direction, and using 
the equations of motion (|2.17|) . All of the non-quadratic dependence is contained in the 
interaction terms Sint — 0{v 3 , z^ 4 ) which are suppressed by powers of g 2 . Here S2 is the part 
of the action quadratic in the fluctuation field, which can be expressed as 



S 2 [u) 



d 3 x v(x) i>(x) 



d 4 x v(x) A [x\ v(x) 



(4.6) 



in terms of the operator A of quadratic fluctuations in the 0° background, 

ss 



5*(x-y)A[x] = - , , 

b<p{x) b<p{y) 

Expanding the boundary functionals gives 



S\x-y){dl + V"W)} . (4.7) 



r f + vf 



i / d 6 k u f {k) 4> c f {-k) + B f h*,(j) c f + B f [0,u f 



after using the boundary condition (|2.21|) , and 

Bi[0,tf+Ui] = 

1 ( d 3 k Ui{k) ti(-k) + Bi[Q,tf] + Bi[0,Ui 



(4.9) 



an ■ v% 

OCR ■ <j>i 



Ot-L ■ Vj 



after using the boundary condition fl2.1S| ). The 4>ij terms in ( |4.8| ) and (O) cancel with the 
boundary terms in the action (|4.5|) . Note the appearance of additional terms linear in in 
( |1.9| ), which arise because we are expanding around a stationary point of the effective action 
( |2.16 ), instead of the action. Note also that the terms independent of Uij are the overlaps 
of the classical field with the initial and final states, respectively. 
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Now, substituting flOf) , (fO ' 



) and flOf) in ( ^4|) , we obtain 



5[b*,2] = / d/ , e %[b^] + Bi[oA c ] + ^[^] Z ^ C ] + (41Q) 

where the quadratic integral over fluctuations is contained in 

Z W\ = (4.11) 
J dvfdv.Dv {1 + a R -Vi} {1 + ot^-i/i} e %[o^/] + -Bi[o^]-^-^-^^+^[ I /] ) 

after using • 0? = 1 in the initial state factors and in the exponent. An integration over 
collective coordinates has been factored out of Q4.10] ) and the remaining functional integral 
Q4.ll ) involves only fluctuations orthogonal to any zero modes of the classical background. 



The 0(g 2 ) corrections result from the interactions in S int . Expanding exp^S^fz/]) inside 
( |4.4|) yields the leading term plus corrections which have the form of vacuum to vacuum loops. 
The vertices used to construct these loops each carry a suppression of g 2 . The propagators 
inside these loops are the usual Feynman propagators evaluated in the classical background 
(f) c i. In the instanton case [fj, [TTJ , dangerous initial state corrections arose from the residues 
of propagators in the instanton background, which displayed an energy dependence ~ g 2 s, 
where s is the center of mass energy of the collision and could be parametrically of order l/g 4 . 
Note that it is only the residue of the instanton propagator that is ill-behaved and not the 
propagator itself. Since the corrections in our case are proportional to loop integrals rather 



than the residues of individual propagators 0, [TTJ] , we expect them to remain subleading 
even at high energies. 



The Gaussian path integral appearing in ( [4.1 1| ) can now be done exactly. We first factor 



out the quadratic prefactor by defining a more general functional integral, 

Z [j] = J du f dvi Dp e B f [0 '"f ] + + lS * [u] ' $ d * k , (4.12) 

in terms of which ( fOT) may be written 



ZW\ = jl + j 1 + aL -j]\ Z ^ I , (4-13) 



13=3* 



where we will set the arbitrary current j(k) equal to 

j*(k) = V2^ (a R (k) + a L (k)) e-^ Ti , (4.14) 



6 It is worth remembering that since we are working directly in Minkowski space (instead of beginning in 
Euclidean space and analytically continuing back), the Feynman boundary conditions on Greens functions 
are only obtained through some form of regularization, such as an ie prescription. 
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after the functional differentiation. The stationary point of (|4.12|) is a "classical" fluctuation 
field v c satisfying 

A [x] v c (x) = , (4.15) 

with boundary conditions 

iu^(k) + u k u°(k) = j(k) and -iv c f (k) + LO k v c f (k) = 0. (4.16) 

The solution of this linear homogeneous boundary value problem can be expressed in terms 
of a Greens function 

u(x) = J d 3 qG(x,q;t)j(q) , (4.17) 

where G obeys 

A[x] G{x,x'-t) = , (4.18) 
with boundary conditions in terms of its Fourier transform 

2G(£,g;X t ) + uj k G(k,q;T l ) = 5 3 (k - q) (4.19) 

and 

-iG(k,q;T f ) + uj k G(k,q;T f ) = 0. (4.20) 

Using the above expressions, we obtain a compact expression for ( [4.12 ) entirely in terms 
of this Greens function, evaluated at t — T i; and the known initial state source. 

Z [j] = det- 1/2 [A] exp Q J d 3 k d\ j(k) G (k, q; T t ) j(q) ) . (4.21) 

The exponent appearing here is due entirely to the boundary functional Bi in ( [4.12 ). All 
other terms vanish on the stationary phase solution ( 4.17| ). With this result, we obtain the 
Gaussian integral ( 4.13Q by functional differentiation and setting j = j* 

Z[<j) c } = det~ 1/2 [A] . (4.22) 

Several terms in ( 4.13| ) have vanished in the limit of Tj — > oo, most notably initial state 
factors. This is because both j*(k) and our inner product ( |2.14| ) contain rapidly oscillating 
factors e~ tuJkTi . This rapid oscillation allows the application of the Riemann-Lebesgue lemma 
in the limit T{ — > — oo, which guarantees that for any function f(k) whose Fourier transform 
exists, 

lim / dk f(k) e~ iWkT > = 0. (4.23) 

Here the functional determinant det _1,/2 [A] is actually the product of three determinants: 
the standard 4-dimensional determinant of the operator A satisfying Feynman boundary 
conditions, and two 3-dimensional determinants representing the edge fluctuations v(x, Tij). 
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The latter determinants can be expressed in terms of a homogeneous Greens function similar 
to G described above. 



We have demonstrated that the kernel of the S-matrix ( OJ) is obtained to leading order 
in g 2 from the action of the classical solution and the determinant of quadratic fluctuations 
around it. 

S\b*,2] = Jdfi e B &*'*°f ] + m[0 ^ ] + W] det" 1/2 [A] [l + 0(g 2 )] . (4.24) 

The classical factors here are known functions, depending only on the initial and final bound- 
ary values of the classical field, ( |2.19| ) and ( |2.22j) . Thus, we may further reduce ( |4.24| ) to 

S[b*,2] = e' 1 J dfie N/2 + iS[ * c] det" 1/2 [A] [l + 0(g 2 )] . (4.25) 

where N = J d 3 k b^b^ is the average number of particles in the final coherent state, which 
depends implicitly on the classical field. The origin of this factor is the use of coherent states 
which are not normalized. The standard normalization for coherent states which we use here 
is 

(b*|b*) = e^. (4.26) 

The correctly normalized S-matrix element would not include this factor, leaving the classical 
action and the determinant. 

There are no additional tree-level radiative corrections to ( [4.25] ), as there are in an expan- 



sion around a stationary point which does not obey the correct boundary conditions (c.f.- 
the instanton [0, H]). This is a direct result of the fact that our classical solution obeys both 
the correct equation of motion and the correct boundary conditions. 

It is worth noting that our semiclassical expansion differs from the usual one in a crucial 
way. Our classical background depends on the coupling constant itself. In the usual instanton 
calculation, the field equations and boundary conditions can be made independent of g by 
rescaling <p — > <f) = g<ft. In our case, since the boundary conditions are fixed by the wave 
packet shapes otR^ik), the equations of motion for contain a dependence on g which cannot 
be scaled away. It is clear that a change in g will lead to a different nonlinear evolution of 
4> and hence a different classical trajectory. In the limit of strictly zero coupling constant, 
the evolution of the initial wave packet is clearly trivial, and we get only a contribution to 
the "diagonal" part of the S-matrix (i.e.- no scattering). Therefore, in order to obtain a 
nontrivial result, we have to take g large enough to find an interesting classical trajectory, 
and yet small enough to justify the semiclassical g 2 expansion derived above. We have not 
demonstrated that the two requirements can be satisfied simultaneously in any theory. This 
can only be checked a posteriori by computing the size of the subleading 0(g 2 ) corrections 
in the background of a proposed trajectory. Of course, our experience with perturbation 
theory leads us to believe that for g 2 <C 1 the corrections will be small. Thus, a classical 
trajectory found using weakly coupled equations of motion is likely to provide a controlled 
approximation. 
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5 Complex Trajectories 



In the preceding discussion, we have focused on classical trajectories dependent on a real time 
variable. These are relevant to classically allowed amplitudes, such as those for multiparticle 
scattering or transitions at energies well above an energy barrier in configuration space. 
As we noted previously, the existence of a real time trajectory implies an unsuppressed 
scattering amplitude for the corresponding transition. However, there are many instances 
in which classically disallowed transitions are of interest. In particular, baryon number 
violation in the electroweak theory at energies less than the sphaleron energy E s ~ M w /a w ~ 
10 TeV is classically forbidden and therefore cannot be described by the same type of real 
time trajectories. The application of our formalism to such cases therefore requires some 
generalization, which we will describe in this section. 

Consider the initial value problem given by ( p,17|) and Q2.18|) . Since both the initial field 
values 4>i(x) and their time derivatives <pi{x) are given, the classical field has a well defined 
energy. Since real time trajectories conserve energy, it is clear that a suitable trajectory will 
not exist for some choices of initial and final conditions. In particular, in the electroweak 
theory for E c i < E s , there are no real time trajectories which connect initial and final 
conditions with different topological charge. 

Im t 



iT 



-> Re t 



Figure 2: Complex time contour 



T2], O, IT1] have proposed complex 



In order to find such trajectories, previous authors 
time paths, which involve alternately Minkowskian, Euclidean and Minkowskian paths joined 
together at two times, say and iT as in figure 2. In this approach, t = and t = iT must be 
turning points (<fr(x, Tj) = for all x) of the field equations with given initial conditions. As 
is familiar from one dimensional quantum mechanics ||13|| , complex time contours describe an 
incident configuration which reaches a turning point, and then tunnels under the barrier in 
Euclidean time, returning to real time at the second turning point upon reaching the other 
side of the barrier. The field d(x) can be taken to be real on the entire contour because 
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of the existence of the two turning points @. The action is then purely imaginary on the 
Euclidean part of the contour and yields the familiar exponential tunneling suppression. 



</IO 



-2lmS 



(5.1) 



As a matter of practice, one might try to find a suitable complex time trajectory by 
integrating our initial conditions (|2.20|) forward in real time in hopes that a turning point 
is encountered. Then, the time variable should be Wick rotated t — > ix^ and the turning 
point configuration integrated forward in Euclidean time until another turning point is en- 
countered. Finally, the contour should be rotated back to real time and the second turning 
point integrated forward until it achieves its asymptotic plane wave state. The final state 
| b*) can be read off from ( |2.22|) . 

However, in the most general case, there is no guarantee of encountering turning points. 
Here we propose a more general method for extremizing S-matrix elements for which there 
is not a corresponding real, Minkowskian trajectory. Rather than complexifying time, we 
instead search for a general complex saddlepoint configuration satisfying ( |2.17|) and ( 2.18|) . 
Our motivation stems from the usual method of steepest descents applied to a real integral 
of the form 

1(h) = [ " dx el f{x) . (5.2) 



Im x 




In general, the asymptotic series in h for this integral can be obtained by expanding 
about a saddlepoint of the function f(x) : f'(x)\ x * = 0. However, the saddlepoint x* is often 
complex, and to apply the method, one must first deform the original real integral into the 
complex plane (see figure 3). In our case, we first think of the Feynman Path Integral as the 

A fact guaranteed in one dimensional QM, but not in higher dimensions or in field theory El. 
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product of a large number of regular integrals by discretizing spacetime: 

J D( p e iS[4>] _^ J]] d$ x e «(*-,*«±0. (5.3) 



For a real field <f)(x) each integral in the product is a real integral like ( |5.2|) . We generalize 
the method of steepest descents by allowing the variable of each integration, $ 2 , to become 
complex. There should be no obstruction to this generalization as the function S(<& x , $ x ±i) 
is analytic in the variables <$> x . This procedure leads us to search for a complex trajectory 
(or in the continuum, <p c (x)) satisfying (|2.17| ), ( |2.18| ) and ( |2.21 ). The integral over 



small fluctuations about c (x) must be performed along the path of steepest descent in 
configuration space. In other words, along small v(x) fluctuations which keep the imaginary 
part of iS [4> c (x) + u{x)\ constant. 

For a complex field, there is no longer a relation between the positive and negative energy 
Fourier components. Without loss of generality, we can write 

(p) = ~4= (u p e~^ + v p e^) , (5.4) 
'2lu p 



where u p , v p are independent complex functions. The initial boundary condition fl2.18|) is 
satisfied by 

u p = + M ^ , (5.5) 

/ d 3 k or(&) / d 3 k a>L(k) 



where t>k is arbitrary. The initial configurations which satisfy (|5.4j ), (|5.5| ) have negative 
frequency modes which "resemble" the wave packet ai + ocr (up to complex multiplicative 
factors) and positive frequency modes which are arbitrary. 

It is then clear that there are an infinite number of complex trajectories which can result 



from evolving ( |5.4j ) using ( 2.17|) . The final quantum state |b*) which results from a given 



trajectory depends only on the positive frequency component of the trajectory in the far 



future (see Eq. ( |2.21|) ). The multiplicity of trajectories allows a very large set of scattering 
amplitudes to be addressed within our formalism. In fact, it seems in general there always 
exists a complex trajectory which satisfies ( |2.17| ), (|2.18|) and fl2.21|) simultaneously for any 



a(k) and b*. This is because the initial and final conditions ( 2.18 ) and Q2.21|) each correspond 



to conditions on a combination of 0(x)|^ and 0(x)|^, and 4>{x)\t } and cj)(x)\T f , respectively 
(mixed Dirichlet-Neumann boundary conditions). Since the equations of motion ( p. 17 ) are 



second order in time, this mixed set of boundary conditions should be enough to specify a 
unique solution. 

Since the trajectories found here are complex, the action iS[(j) c } may contain a real sup- 
pression factor. This is to be expected, as some of the trajectories will correspond to clas- 
sically disallowed tunnelling transitions. Thus, although potentially any S-matrix element 
can be approximated in this way, the result for particular trajectories may be exponentially 
small. 
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6 Nostalgia 



Here we comment on the relationship between our approach and some previous work by other 
authors on high energy multi-particle and anomalous processes in electroweak theory. First, 
we will recall the expansion around the constrained instanton in electroweak theory [fj, |J. 
This approach neglects effects of initial and final states on the saddlepoint and results in a 
low energy (E -C E* ~ M w /a w ) approximation to the total two-particle cross-section. Then, 
we remark on a previous approach to account for the impact of initial and final states on 
the saddlepoint || and comment on the formal resemblance which it bears to our approach. 

A few of the properties of the electroweak instanton will make clear the relationship to 
our work: (i) The electroweak instanton is a solution of the Euclidean equations of motion. 
So, it will be necessary to rotate our real time formalism to Euclidean space in the analysis 
above, t — > 2x4. We would then be considering a saddlepoint approximation to an on-shell 
truncated Euclidean Greens function, (ii) The electroweak instanton has finite Euclidean 
action, S = S7i 2 /g 2 , and therefore satisfies vacuum boundary conditions. So, it is a suitable 
saddlepoint point for a Greens function with external fields only insofar as the effect of initial 
and final states on the saddlepoint are neglected (dropping the right hand side of ( |2.18|) and 

O))- 

Since the instanton does not obey the correct boundary conditions to be a saddlepoint of 
a scattering amplitude 0, linear terms in the fluctuation expansion do not cancel, and the 
expansion of the scattering amplitude entails corrections which are formally large 0(1/ g 2 ) in 
the exponent. Fortunately, these corrections are under control and calculable for sufficiently 
small energies, E <C E* M. The effect of the final states can be taken into account, in a 
perturbative expansion in powers of x = E/E* -C 1 0. For instance, the total inclusive 
cross section in the one-instanton sector is given by the well-known result 

crUE) ~ exp |H£ (-1 + Cl x^ + c 2 x 2 + O (* 8/3 )) j , (6.1) 

in the limit E / 'E* = fixed and g — > 0. Here ~ implies that only the exponential behavior 
of the cross-section is shown. The first term here is just twice the instanton action, the 
'tHooft suppression for vacuum tunnelling, while the next term indicates the exponential 
growth of the cross-section at low energies PJ. The higher order terms are determined from 
all tree- level corrections to the many soft (E ~ M w ) final state particles, and the first few 
coefficients Cj are known. Thus, provided the energy is sufficiently low E <C E*, tree-level 
corrections to the final state can be described semiclassically. 

It has not been similarly demonstrated that the effect of initial state corrections, involving 
many loops, can be described semiclassically. However, some results indicate that these 
corrections may run counter to naive intuition and the initial two particle state may be 
described semiclassically as well |Tl], 0,0. As we noted in Section 4, initial state corrections 
in the instanton expansion appear as factors of the residue of the propagator in the instanton 
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background. The residue of the propagator has hard high energy behavior ( ~ g 2 s ) where 
s is the center of mass energy Jll|, so that initial state corrections can contribute at the 
semiclassical level when s is of order 1/g 4 . Mueller ]lTJ has shown that they contribute first 
at order 10/3 to (p]). 



By contrast, in our formalism the classical background field is constructed to solve the 
correct multi-particle boundary value problem derived in Section 2. So, there are no addi- 
tional tree- level corrections in the limit E = fixed while g — > 0. All corrections are formally 
suppressed by powers of g, as shown in Section 4. The size of subleading corrections to 
our result depend on the properties of the Feynman propagator in the background of our 
proposed classical trajectory, <ft c . It would therefore be very interesting to understand the 
high energy behavior of the propagator in our background field. Also, for the purposes of 
direct comparison to ( |6.1| ), it would be interesting to obtain its behavior in the same limit, 
where E/E* = fixed while g — > 0. We have reason to expect that the propagator behaves 
less severely than the instanton propagator in this limit, due to the lack of translational zero 
modes of the background field ||. 

A previous approach to account for the impact of initial and final states on the saddle- 
point || bears a resemblance to ours but results in a different classical equation. We discuss 
it here for the purposes of comparison. These authors extremized Minkowskian n-point 
Greens functions of electroweak gauge fields 

G (n) { Pl ,p 2 , ...) = JvA A" iai ( Pl ) ■ ■■A ftna "(p n ) e iS[A]/g2 . (6.2) 

It is not clear to us whether the Greens function is the proper quantity to extremize to yield 
information about the corresponding scattering amplitude (and ultimately the cross-section). 
It may be that extremizing the entire Greens function, rather than its LSZ projection, is too 
strict a requirement. 

An equation of motion for the saddlepoint of ( |6.2|) can be derived in a manner with 
formal similarity to our own. To do so, first express the n fields as an exponential, using 
A = explnA, and then extremize the exponent. This yields || 

■ V 8^ 8 aib A -. . . (6.3) 



5A ub (x) t A^ a i(pi 

This is a non- linear integro-differential field equation, depending on both the field A(x) and 
its Fourier transform A(p). Not much is known about equations of this type, and despite 



much effort (|6.3|) defies solution fl5|| . It appears more difficult to solve in practice than the 
boundary value problem we have derived ( |2.17|) , ( |2.19|) and ( |2.21| ). Of course, we do not yet 
have an explicit solution to our boundary value problem, but we hope to have established a 
practical procedure for computing one. 
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7 Conclusions 



In this paper, we have established a connection between certain classical trajectories and 
corresponding quantum S-matrix elements in a weakly coupled field theory. The connection 
is made via a semiclassical approximation which appears to be controlled at small coupling, 
regardless of the scattering energy. This is in contrast to methods which involve expansions 
about instanton-like trajectories (i.e.- satisfying vacuum boundary conditions), which have 
been shown to break down at nonperturbative energy scales E ~ m/g 2 . Real Minkowskian 
trajectories can be found by the straightforward time integration of a well-defined initial value 
problem which is determined by the inital quantum state. Complex Minkowskian trajectories 
require the implementation of mixed boundary conditions at both initial and final times. 
This procedure, while more difficult in practice, can be guaranteed in principle to describe 
any wave packet to coherent state transition specified by the complex functions a(k),bl. 
The formalism described here can thus be used to compute intrinsically nonperturbative 
scattering amplitudes in a variety of field theories. 

Some examples discussed here are the problem of baryon number violation in the elec- 
troweak theory and multiparticle scattering in gauge and scalar theories. The discovery of 
nontrivial or "interesting" real classical trajectories in these theories can now be related to 
unsuppressed scattering amplitues which are potentially observable experimentally. Alter- 
natively, nontrivial complex trajectories, of which there are an infinite number, allow the 
calculation of a much larger class of scattering amplitudes, including some that are clas- 
sically forbidden and therefore involve tunneling. We hope that our results will stimulate 
future numerical work in this area, particularly the search for nontrivial trajectories. 
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